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The energies of the excited states of the Nucleon, A and O are computed in lattice QCD, using 
two light quarks and one strange quark on anisotropic lattices. The calculation is performed at 
three values of the light quark mass, corresponding to pion masses = 392(4), 438(3) and 521(3) 
MeV. We employ the variational method with a large basis of interpolating operators enabling six 
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in the pattern of states. The need to include operators that couple to the expected multi-hadron 
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I. INTRODUCTION 

The theoretical determination of the spectrum of 
baryon resonances from the fundamental quark and gluon 
degrees of freedom is an important goal for lattice QCD. 
Recently a number of groups have performed lattice com- 
putations for Nf = 2-1-1 QCD with two light-quark fla- 
vors (up and down) and one more massive quark flavor 
(strange) . P]-jS] Different actions have been used and sev- 
eral analyses achieve a pion mass close to the physical 
limit. A recent analysis used a reweighting technique to 
perform calculations at the physical pion mass.^ One 
focus of attention has been the determination of the low- 
est baryon mass for each isospin and strangeness. Good 
agreement has been achieved between different calcula- 
tions. The experimental masses are reproduced typically 
with discrepancies ranging from 1% to 8%. 

Another focus of attention is the excited state spec- 
trum of baryons. Recent works have progressed beyond 
quenched QCD (TUl [H] which omits the effects of light 
quarks from the gauge ensembles. Two light-quark fla- 
vors {Nf — 2) were used in Ref. [T^. In this work we use 
ensembles of gauge configurations developed in Ref. [Tj 
for Nf = 2-1-1 QCD with dynamical light and strange 
quarks. 

There are ongoing experimental programs aimed at de- 
termining the spectra and properties of excited baryons 
at the Thomas Jefferson National Accelerator Facility. 
Reference [? ] provides an overview of some recent ex- 
perimental results. 

Excited baryon states can be quite massive and a small 
lattice spacing in the time direction is best for observing 
their signals, which decrease exponentially with time. If 
the same lattice spacing were used in spatial and time di- 
rections the computational cost would be unnecessarily 



high. The Hadron Spectrum Collaboration has under- 
taken a program to solve QCD using anisotropic lattices 
that have a smaller spacing, at, in the time direction. In 
this work we use 16'^ x 128 lattices with a renormalized 
anisotropy <^ = 3.5 [13], i.e., at = 0^/3.5, where is the 
spatial lattice spacing. A Symanzik-improved gauge and 
a clover- improved Wilson fermion action have been used 
to generate gauge ensembles for Nf = 2-1-1 QCD. [1] 
The three ensembles used in this work have pion masses 
of 392(4) MeV, 438(3) MeV and 521(3) MeV. The effects 
of charm, bottom and top quarks are neglected because 
the baryons under study have only light quarks in their 
valence structures and the loop contributions of the ne- 
glected quarks are suppressed by their large masses. 

For families of particles with given isospin and 
strangeness, spectra are calculated in the six double- 
valued irreducible representations (irreps) of the octa- 
hedral group. There are three irreps for even-parity that 
are labeled with a g subscript (gerade) and three for odd- 
parity that are labeled with a u subscript (ungerade). 
They are: Gig,Hg,G2g,Giu, and G2u- Operators 
that transform according to one of these irreps do not 
mix with those of other irreps because of the octahedral 
symmetry. 

A large basis of interpolating field operators is needed 
in order to extract the spectrum of excited states. We 
have developed many such operators for each isospin, 
strangeness and octahedral irrep in previous works. |T5] 
In this work we select sets of 7 to 11 operators in each 
irrep. Continuum values of total angular momenta show 
up in lattice simulations as patterns of degenerate ener- 
gies in the continuum limit that match the patterns in 
Table |l] for the subduction of spin J to the double- valued 
irreps of the octahedral group. For example, a state in 
one of the G2 irreps is a signal for the subduction of 
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TABLE I: The number of occurrences of double-valued irrep 
A of the octahedral group for half- integer values of continuum 
spin J. N is the dimension of the irrep. 
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continuum spin | or higher. Because spin | has six mag- 
netic substates and the lattice G2 irrep has dimension 
two, there must be four degenerate partner states from 
the four-dimensional H irrep in order to realize the six 
linearly independent states that are required. For spin |, 
there must be degenerate partner states in the Gi, H and 
G2 irreps in order to realize the total of eight magnetic 
substates. 

This paper is organized as follows. In section |llj we 
review the lattices used and the basic features of the ac- 
tion and quark masses, the baryon operators, the smear- 
ing of quark fields based on eigenvectors of the gauge- 
covariant lattice Laplacian, the pruning of operators to 
obtain suitable sets for calculations of matrices of correla- 
tion functions and the variational method with eigenvec- 
tors fixed at one time value. As we have shown in previ- 
ous work |12j , smearing of the source and sink operators 
is important for calculations of baryon spectra because 
it reduces contributions from short wavelength fluctua- 
tions. We use the recently developed distillation method 
of Ref. [16| in which source and sink quark operators are 
smeared by applying eigenvectors of the gauge-covariant 
Laplacian, — V^. 

The orthogonality of the eigenvectors produces a fac- 
torization of the problem. One part consists of the 
smeared baryon operators based on A'oig eigenvectors of 
the scalar Laplacian — on each time slice for each 
quark. The other consists of correlation functions cor- 
responding to the parallel transport of smeared quark 
operators formed from the eigenvectors of the Laplacian 
from the source time slice to the sink time slice. These el- 
ementary quark correlation functions are called "peram- 
bulators" for short. They are universal in the sense that 
they may be used for any baryon or meson operators. 
The distillation method provides an 'all-to-all' calcula- 
tion in the sense that propagators are calculated from all 
iVoig eigenvectors of — at the source to all A'oig eigen- 
vectors at the sink. 

In Section |III| we report calculations for the nucleon 
and show details of the fits to obtain energies and the re- 
sults for the spectrum at the lowest pion mass, 392 MeV. 
In Section [TV] we report summary results for the nucleon, 
A and VL spectra for all three pion masses and compar- 
isons with experimental resonances. The f2 spectrum of 
low-lying excited states is interesting because little ex- 
perimental information is available. It also helps to set 
the overall scale for baryon spectra, for which we use the 



TABLE II: Lattice parameters. Three sets of values (in tem- 
poral lattice units) of the light quark mass, mi, the strange 
quark mass, ttIs and the resulting vr meson, K meson and Sl- 
baryon masses. The corresponding values of in MeV are 
given in the last row based on using Eq. ( |15| |. 



ensemble 


12 3 


mt 


-.0840 -.0830 -.0808 


ms 


-.0743 -.0743 -.0743 


Volume 


16^ X 128 16^ X 128 16^ x 128 


A'cfgs 


344 570 481 


^sources 


4 4 4 


m„ 


0.0691(6) 0.0797(6) 0.0996(6) 


rriK 


0.0970(5) 0.1032(5) 0.1149(6) 


mn 


0.2951(22) 0.3040(8) 0.3200(7) 


(MeV) 


392(4) 438(3) 521(3) 



ground state fl baryon mass at each value of m^. 
Section [V| gives a summary of the results. 

II. LATTICES, OPERATORS AND MATRICES 
OF CORRELATION FUNCTIONS 

References [T] and [13] presented the tuning of quark 
masses and other parameters for Nf^2+1 QCD on the 
anisotropic 16'^ x 128 lattice used in this work. Refer- 
ence |16j presented the method of smearing hadronic op- 
erators using an expansion in terms of eigenvectors of 
the three-dimensional lattice Laplacian. Because those 
references provide the foundation for the present work, 
we summarize their findings in this section. 

Table [IT] shows the three sets of quark masses from 
Ref. [T] used in this work together with values for the 
pion mass, the kaon mass and the baryon mass. The 
value of the pion mass in MeV units is obtained by us- 
ing the mass at each point to set the scale. We also 
show the number of gauge configurations used. For each 
configuration we compute correlation functions from four 
different time sources and average them in order to take 
account of correlations. 

Baryon operators used in this work were developed in 
Ref. [T5]. For the baryons that we have considered, the 
single-site (SS) forms of the operators are given in Ta- 
ble [lllj Most operators incorporate gauge-covariant dis- 
placements of the quarks relative to one another in or- 
der to obtain nontrivial shapes. The displaced operators 
are projected to irreps of the octahedral group by sum- 
ming over all lattice rotations applied to the shapes and 
spins with coefficients that project out the irrep opera- 
tors. Many operators are so obtained and in the final step 
they are "pruned" to sets of between 7 and 11 operators 
for each TV and A irrep. For J7, we use the same opera- 
tors as for A except that the number is limited to 6 in 
each irrep. The operators are selected to have good sig- 
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nal to noise characteristics and correlation matrices with 
low condition numbers at a time close to the source. Low 
condition number (ratio of largest to smallest eigenvalue) 
is a proxy for linear independence, which follows because 
two linearly dependent operators would give a zero eigen- 
value for a matrix of correlation functions and an infinite 
condition number. A full listing of operators is available 
upon request. 



TABLE III: Baryons and the corresponding three-quark ele- 
mental operators. Columns 1 to 4 show the symbol, isospin, 
strangeness, and the form of elemental single-site operators 
used. The last three columns show the numbers of embed- 
dings of single-site operators with irreps G\g,G2g and Hg, 
respectively. 
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Table IV indicates the number and type of pruned op- 
erators for each irrep for N , A and baryons used in 
this work following the conventions of Ref. [TBi . The 
types of operators are: SS (single-site); SD (singly dis- 
placed) with one quark gauge-covariantly displaced from 
the other two; DDI (doubly-displaced-I) with two quarks 
displaced in opposite directions from the third; and TDT 
(triply-displaced-T) with all quarks displaced to create a 
T shape. 

In order to reduce couplings to short-wavelength lat- 
tice fluctuations, smearing of the operators is performed 
as in Ref. [TB]. The method is called "distillation" and 
it uses the eigenvectors of the gauge-covariant, three- 
dimensional Laplacian operator. The k*'' eigenvector de- 
pends on a color index, a, and the spatial coordinates, x 



and is written as Va 



(fc) 



It obeys the eigenvalue equation 

(1) 



2\<^b ^^(k) _ . (k) 
b.y ^ '^k^avL 



where a and h are color labels. Ordering is imposed such 
that increasing k corresponds to increasing eigenvalues 

A sum over all M = N^. x N,j. x Ny x eigenvectors 
for a given lattice provides a decomposition of unity on 
time slice t, i.e., 



M 



(2) 



A sum over only the A^eig lowest eigenvalues provides the 
distillation operator on time-slice t as the following pro- 
jection. 



(3) 



TABLE IV: Numbers of operators of each type used in this 
work for N, A and Q matrices are listed for each irrep. 
SS denotes single-site (local) operators, SD denotes singly- 
displaced operators, DDI denotes doubly-displaced-I opera- 
tors, DDL denotes doubly-displaced-L operators and TDT 
denotes triply-displaced-T operators. Gauge-covariant dis- 
placements are used. 



Baryon 


Operator 
type 


trig 




Hg 








AT 
iV 


DO 


o 
Z 


-1 
i 


1 
i 


1 


U 


U 


AT 




1 


1 
i 


o 
z 


1 


1 
i 



Z 


AT 


riFiT 


u 


Z 


1 
1 


O 

o 


Q 


z 


AT 
iV 


riFiT 

ULJLi 


o 

z 


z 


A 

4 


Q 
O 


Q 


z 


AT 


1 U L 


z 


A 

4 


Q 

o 


1 
i 


4 


z 


AT 


total 


i 


1 n 
iU 


1 1 


o 


1 1 
ii 


Q 




/\ 


oc 
oo 








I 


2 


n 



u 


A 


SD 


3 


3 


2 





4 


2 


A 


DDI 


2 


2 


1 


2 








A 


DDL 


3 


3 


2 


2 


4 


4 


A 


TDT 


2 


1 


3 


4 


3 


3 


A 


total 


10 


9 


9 


10 


11 


9 


n 


SS 








1 


2 








Q. 


SD 


3 


3 


2 





4 


2 


n 


DDI 


2 


2 


1 


2 








n 


DDL 


1 


1 


2 


2 


2 


4 


n 


total 


6 


6 


6 


6 


6 


6 



Applying the projection to a quark field yields 



where the smeared field operator is defined by 



(4) 



(5) 



The smeared field operator involves a sum over repeated 
indices h for color and y for space. It depends only on 
the Dirac index, a, the Laplacian eigenvector label, k, 
and time, t. It is of rank Aoig x N^, vector on each time 
slice, where N„ is the number of spinor components, i.e., 
four. A similar projection is used for displaced operators 
as discussed in Ref. 16 . 

Matrices of correlation functions are calculated as fol- 
lows. 



c,,(t,i') = E(^*(^'*)^](y'^'))' 



(6) 



fe=i 



where, using single-site operators, 

i?.(x, t) = Cf-'e^'^q'^f- (x, t)q'/^ (x, i)?^^^ (x, t). (7) 

Superscripts a, h and c are color indices while /i, /2 and 
are quark fiavor indices. Constants (7"'^''' weight the 
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various Dirac components as required to form an irrep of 
the octahedral group. When each quark field is projected 
to the space of the iVeig lowest eigenvalues of the covariant 
Laplacian as in Eq. El), the correlation function becomes 
(with sums over repeated indices understood) 



where 



(8) 



(9) 



averaged over configurations to obtain the eigenvectors 
as follows, 

C\,it\0)ul"\t*) = X.Jt*)C,,{to,0)u^;\t*), (11) 

where to is the normalization time for which the eigen- 
values obey A„(io) = 1- Matrix indices i and j each take 
Nop values where iVop is the number of operators used. 
For each gauge configuration, we then calculate matrix 
elements of Cij{t, 0) in the basis of fixed eigenvectors de- 
termined at time t* , defining the ensemble of effective 
eigenvalues as follows, 



Kit) 



(12) 



Further reduction of the correlation functions involves 
contractions of smeared quark fields that have the same 
quark flavor. That is different for each baryon but in- 
volves the same set of contractions as for unsmeared 
fields. With the smearing used here, each nonvanishing 
contraction yields a "perambulator" in the terminology 
of Ref. Hg, i.e.. 



kk 



{t,t') 



- {^^\t)t\t'))A 



.(fc)t 



yz.cxa. 



{t,t')vm'), (10) 



where M is the Dirac matrix and {■ ■ ■)a denotes eval- 
uation with a single gauge configuration. The peram- 
bulators are matrices in A'oig x dimensions for each 
pair of source and sink times, t' and t. As shown in 
Ref. [IS], iVeig — 32 provides smearing comparable to 
that based on Gaussian smearing for 16^ x 128 lattices. 
We use A'oig — 32 in this work. The perambulator matri- 
ces provide the quark propagation from all eigenvectors 
of the Laplacian at the source time to all eigenvectors at 
the sink time without reference to the operators that are 
used. 

Matrices of correlation functions Cij{t, t') are obtained 
for four time slices t' . We translate each of these corre- 
lation functions to t' = 0, yielding Cij(t,0), and aver- 
age them for each configuration. The average over gauge 
configurations of Eq. ^ is calculated using the baryon 



operators ^fkimi^) ^'^'^ ^'jJclfh^^^ together with peram- 
bulators T^a(i,0) for the relevant contractions. 

The variational method [T71 [THj is used to extract the 
energy levels. We have exploited the fact that the matri- 
ces of correlation functions are real-valued within noise 
after removing a time-independent phase from each op- 
erator _Bj(x, t) and the corresponding complex-conjugate 
phase from Bj{y,t). Imaginary parts of the matrix el- 
ements are therefore dropped and that helps to reduce 
the overall noise. The hermitian matrices Cij{t,0) be- 
come real-symmetric matrices. 

The generalized eigenvalue problem is solved at time 
t* to obtain eigenvectors of the correlation matrices. Be- 
cause matrix elements of the Hamiltonian involve an av- 
erage over configurations, we use the correlator matrix 



At time t — t* , the average over configurations of A„ (t) is 
equal to A„ ( t) b ecause the average correlator matrix was 
used in Eq. (11). At time t = to both the average A„(t) 



and A„ (t) equal one for the same reason. The use of fixed 
eigenvectors provides a smooth time dependence in the 
diagonal correlation functions, A„(t), and is consistent 
with the fact that the eigenvectors of the Hamiltonian 
should be independent of time. Differences between the 
fixed and exact eigenvectors contribute at second order 
to the difference between A„(t) and A„(t) owing to the 
variational nature of the calculation. The accuracy of the 
fixed eigenvalues analysis has been checked by compar- 
ing with the exact eigenvalues. The two analyses agree 
within uncertainties. 

There are two other reasons for using fixed eigenvec- 
tors. One is that the eigenvalue problem generally be- 
comes ill-conditioned at late times owing to the expo- 
nential decrease of A„(t) gc e~^" as t becomes large. 
The higher energy states tend to zero fastest so the con- 
dition number of the A'op x A'op matrix increases expo- 
nentially cx e(^'''K''~^i°"^(*~*''\ where i^high and i?iow are 
the largest and smallest energies obtained from A^op op- 
erators. For the Aop = 7 to 11 operators used, the con- 
dition number becomes very large by i « 15. Backward 
propagating baryon states can become significant also at 
large times. When that happens the smallest eigenvalue 
can pass through zero and become negative. The use of 
eigenvectors fixed at a time t* significantly earlier that 
i « 15 avoids the conditioning problem. 

The second reason for using fixed eigenvectors is that 
at early times the eigenvectors generally are contam- 
inated by contributions from states above the energy 
range that is determined by Aop operators. By choos- 
ing time ^0 to be as large as possible, one reduces the 
contamination from higher states in the spectrum, which 
are suppressed by factors involving e~^*° as shown in 
Ref. [H] . We usually set t* = io + 1 in order to diagonal- 
ize a matrix that is effectively Ri 1 — H + ■ ■ ■ , with 
eigenvectors dominated by the Hamiltonian, H. More- 
over, a long fitting interval [U^tf) is required for accurate 
determination of the energies. We balance the require- 
ments by choosing to be large while keeping t* small 
enough to allow a well-conditioned determination of Aop 
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eigenvectors. The effective eigenvalues A„(i) are then fit 
as described in the next section in order to extract ener- 
gies. 

III. ANALYSIS OF NUCLEON SPECTRA AT 
= 392(4) MEV 



Spectra have been calculated for the isospin ^ 



(N) 



states, isospin | (A) states and strangeness —3 (Q) states 
using lattices with the three sets of quark masses shown 
in Table [ll] In this section we present a detailed discus- 
sion of our analysis for the N* states at the lightest pion 
mass to show the procedure we employ and to illustrate 
the quality of the data; the procedure for other cases is 
similar. 

Table |V] shows the results of fitting each nucleon di- 
agonal correlation function, A„(i), by a two-exponential 
function of time as follows, 



XMt) = {l-A)e 



-E{t-ta) 



+ Ae 



-E'{t-to) 



(13) 



where E < E' so that E is the energy of interest at large 
time. The second exponential term serves to model the 
contaminations arising from higher-energy states at early 
times. The fit window (ti,tf) is chosen such that coeffi- 
cient A is small. The values of ti and tf used are given 
in Table |V] as are the values of to and diagonalization 
time t* . Although the lowest five energies are shown, 
the number of operators used is 7 to 11 in each irrep as 
shown in Table |IV| Fits are performed for a jackknife 
ensemble of diagonal correlation functions calculated as 
in Eq. (12), producing a jackknife ensemble of fit ener- 



gies whose mean and standard deviation are given in the 
table. 

Plots of the nucleon effective energies, calculated as 



2 lA(t + l) 



(14) 



in Figure [T] for the Gig and Gi„ irreps. Fig. [2] 
g and Hu irreps, and Fig. [3] for the G-zg ana 



are shown 
for the Hg 

G^u irreps. These plots show the values of £^eff obtained 



from Eq. ( 14 ) as vertical bars and EcG calculated using 
the fit function of Eq. (131 in place of \{t) in Eq. (14) 



as dashed lines. Comparison of the dashed lines with the 
bars from the lattice ensembles shows the usefulness of 
two-exponential fits. The term Ae~^ (t-to) models the 
contributions of higher energy states at early times al- 
lowing the exponential term (1 — A)e~^'*~*°' to be de- 
termined over a larger fit window {ti,tf) than would be 
possible using a single exponential. Fit energy E and 
uncertainty of the fit energy, a, are shown by dashed 
horizontal lines ai E + a and E — a extending over the 
fit window. Note that the fits over a long time interval 
provide smaller uncertainties when compared with the 
variations of the effective masses. The latter have contri- 
butions from higher states at early times and local fiuc- 
tuations at late times because they are calculated from 



TABLE V: 
MeV. 
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the correlation function at next-to-nearest times. Note 
also that the statistics allow credible determinations of 
six energy levels in each irrep. 

The same process has been used to obtain N, A and 
n energies at three values of m^. The results are given 
in summary form in the next section. 
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FIG. 1: Nucleoli Gig effective energies are shown for the lowest states in the upper six graphs. The effective energy increases 
from left to right along the first row and continues to increase from left to right along the second row. The lower six graphs 
show nucleon Giu effective energies increasing in the same pattern. Calculations are for m,r = 392(4) MeV. Vertical bars show 
the effective energy and the curved dashed line shows the effective energy calculated from the fit function. Horizontal dashed 
lines show the fit results for E iL a and their extent shows the fitting interval {ti,tf). 
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TABLE VI: Multiparticle thresholds on the lattice for isospin 
/ and strangeness S are shown for each value of in MeV. 
The thresholds are based on the sum of energies of the parti- 
cles with no interactions. 
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IV. N, A AND O SPECTRA AT THREE PION 
MASSES 



The goal of determining the spectra of baryon reso- 
nances from lattice QCD requires an increasing elabo- 
rate analysis as the limit of physical pion mass and large 
volume is approached. Although all lattice states have 
discrete energies at any finite volume, the energies corre- 
spond to single-particle states, interacting multi-particle 
states and mixtures thereof. At a minimum, one needs 
to resolve all the states up to some energy and iden- 
tify them as predominantly resonances or predominantly 
scattering states. The repulsion or attraction of multi- 
particle energy levels at finite volume can be related 
to the momentum-dependent phase shifts; the (model- 
dependent) resonance parameters can then be extracted 
through, say, a Breit-Wigner fit to the phase shift. It 
also should be noted that experimental resonances gen- 
erally involve mixtures of single-particle states and multi- 
particle states and in some cases there may be a linear 
combination of multi-particle states that produces fea- 
tures similar to those of a resonance. 

Although we have spectra for three values of m^r, we 
cannot clearly delineate multi-particle states in the spec- 
trum and are unable to obtain the energy-dependent 
phase shift; that analysis must await the introduction of 
a broader basis of operators. In the following, we do not 
attempt to perform a chiral extrapolation on the spec- 
trum. The couplings of the excited states are in general 
unknown, and we are performing calculations in a region 
in which, as we will see below, multi-particle contribu- 
tions are expected. 



A. Nucleoli spectra 

Spectra for isospin | states (N states) are summarized 
in Figure |4] for each lattice irrep and for = 392, 438 
and 521 MeV. We also show the isospin ^ two-star, three- 
star and four-star experimental resonances with val- 
ues that have a subduction to the lattice irrep. Experi- 
mental resonances [20] are shown by boxes in columns la- 
beled by their values with the height of the box equal 
to the full decay width of the resonance. An inner box 
(color aqua) shows the uncertainty in the Breit-Wigner 
resonance energy. The lattice results are shown as col- 
ored boxes with height equal to 2cr in the columns labeled 
by values of , where a is the statistical uncertainty of 
the fit energy. The lattice energies have been converted 
to MeV units by the formula, 



E = 1672.45 



Eat 
mnat 



latt 



(15) 



Here the ratio of a lattice energy and the lattice mn is 
calculated for each value of and then is scaled by 
the empirical mass of the Jl-baryon, 1672.45 MeV. Tri- 
angles to the right of the lattice spectra for each value of 
in Fig. |4] point to the threshold for scattering states 
(multi-particle states) at that value of ttt-tt . Most of the 
excited states have energies higher than the thresholds 
for scattering states. 

As an example, we discuss the Gig plot, shown in 
the upper left hand panel of Figure |4] The experimen- 
tal spectrum contains three low- lying states together 

7 + 

with one ^ state. For the lattice calculation, we show 
the six lowest energy states at each pion mass. Each ex- 
perimental state shown has a subduction to isospin i. 
Gig, and should occur in the lattice QCD spectrum for 
each value of to^. The lattice states should correspond 
to subductions of continuum states. The relevant contin- 
uum states consist of the experimental resonances and 
scattering states with values that have subductions 
to the lattice irrep. Thresholds for scattering states are 
shown in Fig |4] and are listed in Table |Vlj The listed mul- 
tiparticle states are assigned to lattice irreps following 
Ref. [H] . Because the spatial lattice is a cube measuring 
about 1.96 fm on a side, the momentum is restricted to 
discrete values with the smallest nonzero one being 630 
MeV. Consequently, the scattering states with nonzero 
momenta occur at higher energies. 

The general pattern seen in Figure [3] is that lattice 
states have high energies that decrease toward the ex- 
perimental resonances as is decreased. The density 
of lattice states increases with increasing energy, but we 
restrict the analysis to the lowest six states. That we are 
able to extract six energies in each lattice irrep is testa- 
ment to the effectiveness of the smearing procedure based 
on eigenvectors of the lattice Laplacian. Note, however, 
that we cannot distinguish clearly between single- and 
multi-particle contributions to the spectrum in this cal- 
culation. 
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FIG. 2: Nucleon Hg effective energies are shown for the lowest states in the upper six graphs and nucleon effective energies 
are shown in the lower six graphs with effective energies increasing in the same pattern as in Fig. [l] Calculations are for 
rriTr = 392(4) MeV. Vertical bars show the effective energy and the curved dashed line shows the effective energy calculated 
from the fit function. Horizontal dashed hues show the fit results for E ±g and their extent shows the fitting interval {ti,tf). 
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The nucleon ground state shows up as the lowest state 
in the Gig lattice irrep. Its energy decreases in a regular 

manner toward the experimental value shown in the 
column as decreases. The first excited state decreases 
toward the Roper resonance, iV*(1440, I"*") but remains 
above about 1900 MeV for the pion masses used in this 
work. 

In the negative-parity Giu spectra, there are two low- 
lying states at each value of and they tend toward the 
experimental resonances A''(1535, ^ ) and A^(1650, ^ ) 
as ttItt decreases. A number of higher states also tend 
toward the energy of the A^(22G0, | ) resonance or to 
scattering states such {NTr)s-wave, {NT:)p-wave and so 
on. 

Before we discuss the H and G2 states, it is worth 
noting that isolated G2 states do not correspond to any 
physical state. Because G2 has minimum spin |, there 
must be at least six linearly independent components in 
the continuum limit. Each G2 state must have a part- 
ner H state with the same parity in order to have an 
interpretation as a physical state. However, on the lat- 
tice discretization effects can cause the H and G2 partner 
states to have different energies at 0{a?'). 

In the Hg spectra there are five experimental res- 
onances: iV(1720, 1"*") and iV(1900, |^), iV(1680, |^), 

iV(20G0, 1^) and iV(1990, |^). The lattice states tend 
as a group toward these energies as to^ decreases. 

In the iJ„ spectra there are four low-lying lattice states 
near 1800 MeV. The threshold for scattering states is 
near the same energy as this group of lattice states. 
Three low-lying experimental resonances are present: 
N(153G, I"), N(1650, |") and N(1675, |"). 

In our Nf=2 analysis of Ref. [12], we obtained three 
low-lying iJ„ states with larger uncertainties. Otherwise 
the low-lying lattice states agree reasonably well. A test 
of stability was performed by omitting Hu operators to 
obtain sets of iVop = 6, 7, 8 and 9. Spectra were cal- 
culated for each of these and the results show three 
states near 2000 MeV when we use 6 or 7 operators and 
four Hu states when we use 8 or 9 operators. This behav- 
ior suggests that one state is only resolved with the larger 
number of operators. The operators that are responsible 
for the appearance of the fourth state are of the triply- 
displaced-T type. We also have studied the stability of 
the spectrum at = 561 MeV by varying to, keeping 
t* = to + 1 and using all operators. For = 2, 3, 4, 5, 6, 
7 and 8 we observe four low-lying states in the nucleon 
Hu spectrum. The presence of a fourth state is robust 
when all operators are used. 

Four low-lying Hu states are consistent with experi- 
ment if one is a scattering state. We have not previously 
found any evidence for a scattering state with our three- 
quark operators, but they should be present. Work that 
is in progress aims to identify the scattering states by 
using operators designed to couple to them directly. 

The pattern that is used to identify a spin | state 
on the lattice is a pair of states in the G2U and Hu ir- 



reps that become degenerate in the continuum limit. A 
candidate for this pattern is present in our iJ„ and 6*211 
spectra: the lowest G2U state and one of the four Hu 
states at essentially the same energy. The pattern for 
spin-^ is a triplet of Gi, H and G2 states at essentially 
the same energy. There are candidates for this pattern 
in the positive-parity spectra. However, the presence of 
scattering states makes a secure identification difficult. 



B. A spectra 

The A spectra are shown in Fig. [sj The general fea- 
tures are the same as for the nucleon spectra: the lattice 
states are high and they tend toward the experimental 
resonances as decreases. The A(1232, ground 
state appears as the lowest state in the Hg spectra. The 
next higher Hg state is close to the ground state but 
appears to tend toward the A(1600, |^) state in the ex- 
perimental spectrum. 

The lowest two Gig states near 2200 MeV appear to be 
somewhat high but consistent with the experimental res- 
onances A(1910, i^) and A(1950, |^). Candidates for 
7 -i" 

spin-^ partner states are present in the 2200 MeV to 
2400 MeV range in Hg and G2£, spectra but the pattern 
is a weak match for the expected degeneracy in the con- 
tinuum limit. Possibly the small volume used is causing 
large splittings. 

One of the two lowest iJ„ states corresponds reason- 
ably well to the A(1700, | ) resonance. The other one 
should correspond to the A(1930, | ) resonance. How- 
ever, a suitable partner state for spin | is not seen in 
the G2U spectrum: the lowest such state is near 2600 
MeV. A similar result is found in the G2g spectrum with 
the lowest state being close to 2300 MeV, well above the 
energy of the A(1905, 1"*") resonance. This suggests that 
the volume may be small, particularly for the G2 lattice 
states. In quark models [22l[23], excited states typically 
have larger radii. Our lattice is about 1.8 fm in extent 
and a state with a radius of more than Ifm is not ex- 
pected to be determined well. 



C. fl spectra 

The spectra for excited ft states are shown in Fig. |6] 
The il(1672, |^) ground state has been used to set the 
scale for baryon masses so is reproduced perfectly. Exper- 
imental resonances above the ground state do not have 
spin-parity assignments. In the quality rating of reso- 
nances of Ref. j20j, il(2250) is rated as a three-star res- 
onance while $7(2380) and 17(2470) are rated as two-star 
resonances. The strange-quark mass is at its physical 
value in our calculations and the dependence on the pion 
mass is expected to be smaller than for other resonances. 
Consistent with this the overall pattern of excited states 
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FIG. 3: Nucleon G2g effective energies are shown for the lowest states in the upper six graphs and nucleon G2u effective 
energies are shown in the lower six graphs with effective energies increasing in the same pattern as in Fig. [T] Calculations are 
for — 392(4) MeV. Vertical bars show the effective energy and the curved dashed line shows the effective energy calculated 
from the fit function. Horizontal dashed lines show the fit results for E ±a and their extent shows the fitting interval {ti,tf). 
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varies little with . A noteworthy exception is the first 
Hg excited state, whose energy increases from about 1800 
MeV to 2100 MeV as the pion mass decreases. 

We have considered whether the lattice ft spectra can 
provide a useful guide for assignment of spins and par- 
ities. We find 11 strangeness -3 states with energies 
near or below 2500 MeV. Some of those states may be 
candidates for scattering states rather than resonances. 
Thresholds for scattering states are shown for each value 
of ttItt in Figure |6] However we cannot determine whether 
or not our spectra contain scattering states. 

A reasonably good agreement between the lattice and 
experimental spectra is obtained if the first excited exper- 
imental resonance is assigned to | . Beyond that there 
are several possibilities. In Fig. |6] all of the experimental 
resonances have been shown in the | column that ap- 
pears in the plot of Hg spectra. However, a convincing 
assignment is not possible because many features of the 
spectra are not explained. This issue will be revisited 
when good operators for scattering states are available. 



We find candidates for scattering states that have not 
shown up in our previous analyses based on Nf = 2 QCD 
or quenched QCD. We also find more excited state con- 
tamination in the effective-mass plots than was the case 
for quenched QCD. We expect that when appropriate 
operators are used to identify multiparticle states, the 
spectra will be cleaner. 

The spectrum of excited states of the il baryons has 
been calculated for the first time. We do not find a close 
enough agreement between lattice and experimental ex- 
cited states to allow a convincing assignment of the un- 
known spins and parities. 

Our main conclusion is that the program to determine 
baryon spectra from lattice QCD is expected to produce 
reasonable explanations of the nucleon, A and Q spectra 
once calculations are extended to smaller pion masses, 
larger volumes and operators designed to couple to scat- 
tering states directly. Stochastic estimation of the quark 
propagators will allow use of the distillation method with 
larger volumes. [24] 



V. SUMMARY 

This work represents a milestone in our long-term 
research program aimed at determining the spectra of 
baryons in QCD. It provides the first spectrum for N, 
A and baryons based on Nf = 2-1-1 QCD with high 
statistics. A large number of baryon operators is used to 
calculate matrices of correlation functions. They are an- 
alyzed using the variational method with fixed eigenvec- 
tors. The analysis provides three spectra at pion masses, 
= 392(4) MeV, 438(3) MeV and 521(3) MeV. 

The lattice volume and pion masses used give consid- 
erably higher energies than the experimental resonance 
energies. However, there is reasonable agreement of the 
overall pattern of lattice and experimental states. One 
exception is that almost all G2 states are much too high. 
That may be caused by a volume that is too small for 
highly excited states. 
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FIG. 4: Spectra for isospin | (nucleon family) at three values of are compared with experimental spectra. Plots in the 
first row show Gig and Giu lattice irreps, plots in the second row show Hg and Hu irreps and plots in the third row show 
G2g and G2u irreps. Columns labeled by = 392, 438 and 521 show lattice spectra at those values of m,r. Two, three and 
four-star experimental resonances are shown in columns labeled by their values. Each value listed has a subduction to 
the lattice irrep shown. Each box for an experimental resonance has height equal to the full decay width and an inner box 
(color aqua) showing the uncertainty in the Breit-Wigner energy. Triangles to the right of lattice spectra point to the threshold 
for scattering states at that value of m^. 
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FIG. 5: Spectra for isospin | (A family) at three values of are compared with experimental spectra. Plots in the first row 
show Gig and Giu lattice irreps, plots in the second row show Hg and irreps and plots in the third row show G2g and 
G2u irreps. Columns labeled by = 392, 438 and 521 show lattice spectra at those values of m.^. Two, three and four-star 
experimental resonances are shown in columns labeled by their values. Each value listed has a subduction to the lattice 
irrep shown. Each box for an experimental resonance has height equal to the full decay width and an inner box (color aqua) 
showing the uncertainty in the Breit-Wigner energy. Triangles to the right of lattice spectra point to the threshold for scattering 
states at that value of m^. 
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FIG. 6: Spectra for isospin 0, strangeness -3 {Q family) at three values of are compared with experimental spectra. Plots 
in the first row show Gig and Gi„ irreps, plots in the second row show Hg and Hu irreps and plots in the third row show G2g 
and G2u lattice irreps. Columns labeled by — 392, 438 and 521 show lattice spectra at those values of m^. Two, three and 
four-star experimental resonances are shown in columns labeled by their values. Each value listed has a subduction to 
the lattice irrep shown. For the Q, the spins and parities of the experimentally observed states other than the lightest are not 
clearly determined; for comparison, we assign the states to with the heights of the boxes indicating the widths. Triangles 
to the right of lattice spectra point to the threshold for scattering states at that value of m^. 



